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ABSTRACT 

Galaxy-galaxy interactions are predicted to cause gas inflows leading to enhanced nuclear star 
formation. This prediction is borne out observationally, and also supported by the gas-phase 
metallicity dilution in the inner regions of galaxies in close pairs. In this paper we test the 
further prediction that the gas inflows lead to enhanced accretion onto the central supermassive 
black hole, triggering activity in the nucleus. Based on a sample of 1 1,060 Sloan Digital Sky 
Survey galaxies with a close companion {rj, < 80 /i^,,' kpc, AV < 200 km s"'), we classify 
active galactic nuclei (AGN) based either on emission line ratios or on spectral classification 
as a quasar. The AGN fraction in the close pairs sample is compared to a control sample of 
110,600 mass- and redshift-matched control galaxies with no nearby companion. We find a 
clear increase in the AGN fraction in close pairs of galaxies with projected separations < 40 
hj^ kpc by up to a factor of 2.5 relative to the control sample (although the enhancement 
depends on the chosen S/N cut of the sample). The increase in AGN fraction is strongest in 
equal mass galaxy pairings, and weakest in the lower mass component of an unequal mass 
pairing. The increased AGN fraction at small separations is accompanied by an enhancement 
in the number of 'composite' galaxies whose spectra are the result of photoionization by 
both AGN and stars. Our results indicate that AGN activity occurs (at least in some cases) 
well before final coalescence and concurrently with ongoing star formation. Finally, we find a 
marked increase at small projected separations of the fraction of pairs in which both galaxies 
harbour AGN. We demonstrate that the fraction of double AGN exceeds the expected random 
fraction, indicating that some pairs undergo correlated nuclear activity. We discuss some of 
the factors that have led to conflicting results in previous studies of AGN in close pairs. Taken 
together with complimentary studies, we favour an interpretation where interactions trigger 
AGN, but are not the only cause of nuclear activity. 
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1 INTRODUCTION 

There are several lines of evidence that indicate that the build-up 
of stellar mass in galactic bulges is closely linked to the evolution 
of the central black hole (e.g. Haehnelt, Natarajan & Rees 1998; 
Richstone et al. 1998). Gravitational torques have been shown to 
be effective at funneling gas to the centres of galaxies to feed both 
star formation and accretion onto the central black hole (e.g. Haan 
et al. 2009) and there is a remarkably tight relationship between 
the stellar properties of the bulge and the mass of the black hole 
(Magorrian et al 1998; Ferrarese & Merritt 2000; Gebhardt et al. 
2000). The similarities in the shape of the redshift evolution of 
both star formation rates and the frequency of active galactic nuclei 
(AGN) indicate a common physical driver between the two pro- 
cesses. There is considerable observational evidence linking AGN 
activity with star formation, such as bluer colours in AGN hosts 
(Silverman et al. 2009) and correlations between star formation 



activity and AGN power (Cid-Femandes et al. 2001; Choi, Woo 
& Park 2009). In turn, both observations and simulations suggest 
that feedback from the AGN itself may quench the star formation 
(Kauffmann et al. 2003a; Di Matteo et al. 2005; Schawinski et al. 
2007, 2009; Bundy et al 2008; van de Voort et al. 201 1; Kaviraj et 
al. 2011). 

The obvious culprit for building both the central stellar mass 
and feeding the black hole (and hence increasing its activity) is a 
ready supply of gas for fuel and galaxy-galaxy mergers represent 
a natural mechanism for providing this gas in wholesale quanti- 
ties (Kauffmann & Haehnelt 2000; Cattaneo et al. 2005; Di Mat- 
teo et al. 2005;Springel, Di Matteo & Hernquist 2005; Wild et al. 
2007; Hopkins et al. 2008). Many authors have found evidence 
for mergers associated with AGN hosts (e.g. Surace et al. 1998; 
Canalizo & Stockton 2001; Sanchez & Gonzalez-Serrano 2003; 
lahnke et al. 2004; Smimova et al. 2006; Combes et al. 2009; 
Villar-Martin et al. 2010, 2011), leading to an often cited connec- 
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tion between the merger process and AGN triggering. The discov- 
ery of double AGN also lend circumstantial support to a merger 
origin (e.g. McGurk et al. 2011 and references therein). However, 
although some studies claim a connection between nuclear activ- 
ity and the presence of close companions (Dahari 1984; Keel et al. 
1985; Rafanelli et al. 1995; Koss et al. 2010), other studies have 
claimed that there is statistically no difference in the percentage 
of galaxies with companions between active and inactive galaxies 
(Schmitt2001; Coldwell & Lambas 2006; Grogin et al. 2005). Sim- 
ilarly, there are claims that AGN hosts are no more tidally distorted 
than inactive galaxies (Dunlop et al. 2003; Gabor et al. 2009; Cis- 
temas et al. 2011; Kocevski et al. 2011), but counter claims also 
exist (e.g. Koss et al. 2010; Ramos Almeida 2011a,b). Various bi- 
ases may contribute to these discrepant results, including sample 
size, wavelength-dependent dust extinction, image depth, adequate 
definition of a control sample, automated versus visual classifica- 
tion and distinction between Type I and Type 2 AGN (e.g. Darg et 
al. 2009; Dultzin et al. 2010; Koss et al. 2010; Liu et al. 2011a,b; 
Ramos Almeida et al. 201 lb). Nonetheless, the null results have led 
some authors to conclude that there is actually little evidence for a 
connection between merging and AGN triggering (e.g. Grogin et 
al. 2005; Gabor et al. 2009; Cisternas et al. 2011). These conclu- 
sions potentially undermine one of the cornerstones of our mod- 
ern paradigm of galaxy evolution: the triggering of AGN through 
merger events. 

A natural way to test the connection between mergers and 
AGN is to use close pairs of galaxies which can be classified as 
either star-forming or AGN-dominated on the basis of their emis- 
sion line ratios. Close pairs have been unequivocally demonstrated 
to have enhanced star formation rates (Kennicutt et al., 1987; Bar- 
ton, Geller & Kenyon 2000; Lambas et al 2003; Alonso et al. 2004; 
Nikolic, Cullen & Alexander 2004; Woods, Geller & Barton 2006; 
Woods & Geller 2007; Ellison et al. 2008, 2010; Patton et al. 201 1; 
Liu et al. 2011b) and evidence for gas inflows, as indicated by di- 
luted interstellar medium metallicities (Kewley et al. 2006a; Ellison 
et al. 2008; Michel-Dansac et al. 2008; Kewley et al. 2010; Rupke 
et al. 2010; Scudder et al. in preparation). Similar enhancements 
in star formation rates (SFRs) and low gas phase metallicities have 
been reported for other classes of galaxies that are likely to be post- 
mergers, such as luminous infrared and 'lopsided' galaxies (Rupke, 
Veilleux & Baker 2008; Reichard et al. 2009). It is therefore per- 
haps surprising that the results for AGN fractions in close pairs 
are considerably more controversial. For example, Alonso et al. 
(2007), Woods & Geller (2007) and Rogers et al. (2009) all find 
higher AGN fractions in samples of close pairs, relative to a con- 
trol sample of field galaxies. Conversely, Li et al. (2006, 2008) and 
Ellison et al. (2008) and Darg et al (2009) do not find a statistically 
significant difference in AGN fractions. 

One of the issues that hampers the observational association 
between mergers and AGN ignition is the expected time delay be- 
tween clear signs of merging (e.g. well identified galactic compo- 
nents and strong tidal features), the presence of triggered star for- 
mation and, finally, the activity in the nucleus. This complication 
was highlighted by Storchi-Bergmann et al. (2001) in their study of 
Seyfert 2 galaxies. The majority of the Seyferts with close compan- 
ions show signs of recent nuclear star formation, leading Storchi- 
Bergmann et al. to suggest that the pure Seyfert spectrum emerges 
after the main starburst has faded. A similar evolutionary scenario 
was suggested by Haan et al. (2008), whereby low luminosity AGN 
precede Seyferts after a merger, based on high fractions of the for- 
mer, but few of the latter in a sample of galaxies with disturbed HI 
disks. Delays are also implied by the intermediate colours and spec- 



tral characteristics of AGN hosts, indicating at least ~ 100 Myr be- 
tween the decline of star formation and black hole activity (Schaw- 
inski et al. 2009). Indeed, by the time galaxies have migrated from 
their peak star-forming days into the 'green valley' and/or show 
signs of strong AGN, their morphologies do not exhibit any resid- 
ual signs of interactions (Cisternas et al. 201 1 ; Mendez et al. 201 1). 
Delays of a few hundred megayears between the epochs of peak 
star formation and peak black hole accretion are also inferred from 
simulations (Wild, Heckman & Chariot 2010; Hopkins 2011). Al- 
though gas inflow leads to triggered star formation after the first 
pericentric passage, the highest enhancements are seen when the 
galaxies finally merge. If coalescence also represents the peak in 
AGN activity, this may further contribute to the difficulty in con- 
necting AGN with galaxy pairs. 

In this paper, we use a large sample (Section|2j of close galaxy 
pairs to return to the question of AGN fractions in galaxy-galaxy in- 
teractions. One novelty in the current study is the application of a 
variety of AGN classifications (Section |3j which permit different 
contributions from AGN and star formation. The varying strictness 
of these diagnostics allows us to obtain some measure of how dom- 
inant the AGN contribution is. This will be useful if our sample 
includes a significant fraction of galaxies that are just transitioning 
between their star-forming and AGN phases, the so-called compos- 
ite population. Most previous works have used a single sample of 
close pairs (with separations typically < 30 h^^ kpc) and calculated 
a single AGN fraction to be compared with the field (e.g. Alonso et 
al. 2007; Ellison et al. 2008; Darg et al. 2009). However, one of the 
benefits of our large sample size is that we can investigate the AGN 
fraction as a function of projected separation and as a function of 
pair mass ratio (Section|4j. This has only been attempted once be- 
fore by Woods & Geller (2007) with a sample around 1/3 the size 
of ours. Finally, we will investigate whether there is any evidence 
of simultaneous triggering of AGN (Section |5j. In particular we 
will quantify whether the fraction of double AGN pairs is consis- 
tent with the expectation from random occurrence, or whether there 
is a component of correlated AGN pairs. 



2 SAMPLE SELECTION 

The sample of pairs used in this work differs significantly from that 
used in previous papers in this series (Ellison et al. 2008, 2010; 
Patton et al. 2011) in several respects. Previously, the pairs sample 
was constructed from only the SDSS galaxy (specclass=2) sam- 
ple. However, in this work, we are interested in the AGN contri- 
bution, so we have extended our search for pairs to also include 
objects classified as QSOs (specclass=3) which have been classi- 
fied as galaxies from the SDSS imaging. For clarity, we will refer 
to objects with specclass values of either 2 or 3 as 'galaxies'. The 
parent sample of galaxies from which we will select the pairs is 
restricted to the SDSS legacy area (i.e. excluding the Segue foot- 
print) with extinction corrected r-band Petrosian magnitudes in the 
range 14.0 < m, < 17.77. Our previous works have used stel- 
lar masses calculated from SDSS photometry made available by 
the Max Planck Institut fur Astrophysik/Johns Hopkins Univer- 
sity (MPA/JHU) coUaboratioifl However, in Simard et al. (2011), 
we have shown that the SDSS photometry becomes unreliable for 
galaxies with close companions (see also Figure 10 in Patton et al. 
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201 1). Simard et al. (201 1) recomputed the magnitudes for 1.2 mil- 
lion SDSS galaxies using improved sky subtraction and object def- 
inition techniques, combined with morphological decompositions. 
These improvements dramatically improved the quality of photo- 
metric measurements in close pairs (e.g. Figure 1 1 of Simard et al. 
2011). We therefore re-compute stellar masses using the updated 
photometry of Simard et al. (2011). 

The masses are derived using a simple relation between g-r 
colour and mass-to-light ratio (M/L) derived from MPA/JHU data 
catalogues, which we parametrise as a double power-law. M/L mea- 
surements in the MPA/JHU catalogues were derived from fits to the 
SDSS ugriz model photometry to a suite of synthetic SEDs span- 
ning a broad range of star-formation histories (Kauffmann et al. 
2003b; Salim et al. 2007), and are corrected for the influence of 
emission lines on the SDSS photometry. Model SEDs were con- 
structed using GALAXEV (Bruzual & Chariot 2003) with a Chabrier 
(2003) initial mass function, and galaxies were assumed to be well 
described by the combination of an underlying, exponentially de- 
clining star-formation history and stochastically-sampled bursts. 
Simard et al. (201 1) show that their magnitudes are on average 0.05 
to 0. 1 mag brighter that the SDSS model magnitudes and, more im- 
portantly for this work, are less prone to deblending errors at small 
galaxy-galaxy separations (see also Patton et al. 2011). We there- 
fore adopt our updated masses throughout this work, although note 
that our results are qualitatively similar if we use the MPA/JHU 
stellar masses. 

We impose a minimum redshift of 0.01 to avoid non- 
cosmological values and a maximum value of 0.20 above which 
the sample becomes increasingly incomplete. The final criterion for 
inclusion in the parent sample is that the SDSS redshift confidence 
parameter must have a value of Zconf > 0.7, i.e. the confidence level 
of the redshift is greater than 70 per cent (our results are unaffected 
by increasing this cut to higher confidence levels). 

From the parent sample we select pairs and higher order mul- 
tiplesH with a companion within 80 /i^g kpc, AV < 200 km s"' and 
with stellar mass ratios 0. 1 < Mi I Mi < 10. Fibre collisions lead 
to a high incompleteness at separations < 55 arcseconds which bi- 
ases the mass and redshift distribution of close pairs (Ellison et al. 
2008; Patton &. Atfield 2008). We therefore follow Ellison et al. 
(2008) and exclude a random 67.5% of pairs with 6 > 55 arcsec- 
onds to yield a pairs sample with unbiased selection as a function 
of separation. Before the cull is implemented, there is a clear trend 
towards larger masses at small separations. Since there is a strong 
dependence of AGN fraction on mass, this artificially introduces an 
increased AGN fraction at small separations. After the cull, the dis- 
trubution of masses and redshifts with projected separation is flat. 
There are 1 1 ,060 galaxies in the full pairs sample of which 46 are 
classified as QSOs. 

A control sample is constructed by iteratively matching galax- 
ies in both mass and redshift from a pool of galaxies with no close 
companion within 80 h:^^ kpc and 10,000 km s"' . The control sam- 
ple is vital since the AGN fraction depends sensitively on galactic 
properties such as mass, luminosity, colour, morphology and con- 
centration (Kauffmann et al. 2003b; Best et al. 2005; Choi, Woo 
&. Park 2009; Ellison et al. 2008). Matching in mass largely mit- 
igates these differences, as well as dealing with aperture bias (the 
varying covering fraction as a function of redshift). The matching 
procedure finds the best simultaneous match in mass and redshift to 



^ 94 per cent of our sample consists of pairs, approximately 6 per cent are 
triples and less than 1 percent are higher order multiples. 
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Figure 1. The BPT diagram for control galaxies with strong emission lines. 
Black curves indicate the demarcation lines given by KOI (upper line), K03 
(middle line) and S06 (lower line), see equations[T]to[3] Galaxies are colour- 
coded for easy distinction: star-forming (blue), composite (green) and AGN 
(red). QSOs not present in the MPA/JHU line flux catalogues are excluded 
from this Figure. 

each galaxy (i.e. both members of the pair are matched separately) 
and calculates the Kolmogorov-Smirnov (KS) probability that the 
masses and redshifts of the control and pairs are drawn from the 
same distribution. If the KS probability is above 30%, the match- 
ing is repeated up to the A''* best match. We are able to match 10 
controls to each galaxy in a pair before the KS test fails, yielding 
a control sample of 1 10,600 galaxies. We keep track of which con- 
trols are matched to which pair galaxy so that we can include only 
the relevant controls for different subsamples of pairs. 



3 AGN DIAGNOSTICS 
3.1 Background 

Since their early use in distinguishing spectra dominated by AGN 
(Baldwin, Phillips & Terlevich 1981; Veilleux & Osterbrock 1987), 
line ratio diagrams have evolved into a standard tool for categoriz- 
ing the dominant source of ionizing radiation in a galaxy (Dopita 
et al. 2000; Kewley et al. 2001, 2006b; Kaufi'mann et al. 2003a; 
Stasinska et al. 2006; Groves et al. 2006; Cid-Femandes et al. 
2010). The principle behind these diagnostic diagrams is to dis- 
tinguish the excitation mechanism (photoionization by hot stars, 
AGN or shock heating). The various mechanisms result in a differ- 
ent ionization structure in the interstellar gas, and changes in the 
relative sizes of the ionized, partly ionized and neutral media. The 
inner zone is highly ionized, hosting species such as [OIII]. The 
outer zone is the primary location of the lower ionization species 
such as [Oil] and some of the [Nil]. Forming a partially ionized 
zone requires a relatively hard spectrum as the decreasing cross sec- 
tion of hydrogen (and other species) means that only the high en- 
ergy photons penetrate. This means that the partially ionized zone 
is prevalent in the presence of an AGN and shocks, but virtually 
absent with stellar spectra. [Nil], [OI] and [SII] are all formed in 
the partially ionized zone, and can be boosted further by collisional 
excitation. The differences in structure and ionization balance are 
encapsulated in a single ionization parameter, q: the ratio of hydro- 
gen ionizing photon flux per unit area to the local number density 
of hydrogen atoms. In terms of the choice of line ratios used to 
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Figure 2. Top panel: The fraction of AGN in close galaxy pairs (filled 
points) as a function of projected separation. Open points show the AGN 
fraction for the control sample, where the projected separation conesponds 
to the pair galaxy to which it is matched. Black symbols consider only 
galaxies classified as AGN from the S06 BPT diagnostic. Red symbols 
include the specclass=3 (QSO) objects. Lower panel: The AGN excess is 
determined from the ratio of the AGN fraction in the pairs relative to the 
matched control, as a function of projected separation. AGN are up to 2.5 
times more common in close galaxy pairs than the control (independently 
of whether QSOs are included or not). 



distinguish the ionizing source, preference is usually given to those 
that are relatively close in wavelength in order to avoid a significant 
dependence on internal extinction properties. 

The most commonly used combination of line ratios is 
[NII]/Hq' versus [OIII]/H/3, often simply referred to as the BPT dia- 
gram after Baldwin, Phillips & Terlevich (1981). Although galaxies 
dominated by star-formation and AGN broadly distinguish them- 
selves as two 'wings' on the BPT diagram, the precise demarcation 
is not clear. Using photoionization models, Kewley et al. (2001, 
KOI) proposed a classification scheme that extends into the AGN 
wing of the BPT diagram. The Kewley et al. (2001) classification 
includes the most extreme models of starburst galaxies in terms of 
metallicity and star formation rate. However, galaxies classified as 
star-forming by the Kewley (2001) system may contain up to ~ 
20% contribution from AGN (Stasinska et al. 2006). Stasinska et 
al. (2006, S06) also model a range of ionization parameters and 
metallicities, but with updated model atmospheres, resulting in a 
more stringent demarcation between AGN and star-forming galax- 
ies. Between these extremes, KaufFmann et al. (2003a, K03) have 
proposed an empirical distinction between AGN and Hll-region 
dominated galaxies, based on a sample of ~ 23,000 galaxies from 
the SDSS. 
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Figure 3. The AGN excess (fraction of AGN in the pairs relative to the 
fraction of AGN in the control sample) is plotted as a function of projected 
separation. Black open symbols show galaxies at z < 0.1 and red filled 
points show galaxies at z > 0. 1 . 



3.2 Application to the SDSS sample 

Galaxies in our sample are classified as AGN based on the ratios 
of strong emission lines, or based on their spectral classification 
as QSOs (specclass=3). Emission line fluxes are taken from the 
MPA/JHU catalogues. The catalogue fluxes are already corrected 
for underlying stellar absorption and Galactic extinction. Galaxies 
not present in the MPA/JHU catalogues are discarded (this includes 
some specclass=3 QSOs, most notably those with large Balmer line 
widths). Fluxes are further corrected for internal extinction using 
the ratio of Hq'/^8 and an SMC extinction curve (Pei 1992). The 
minimum criterion for AGN classification is the detection at > 5cr 
of H/3 A4861, [OIII] A 5007, Ho- ^6563 and [Nil] ^6584. There are 
also cuts on the flux error and continuum error of the emission lines 
that clip spurious measurements, see Scudder et al. (in preparation) 
for full details. 

In this paper, we will investigate all three of the common 
BPT classification schemes: Kewley et al (2001), KaufFmann et al. 
(2003a) and Stasinska et al. (2006). The distinction between AGN 
and star-forming galaxies for the 3 schemes are given in equations 
[I](K01),|2](K03) and[3](S06). 

log([OIII]^5007/H/3) = , — T-r; + 1-19 (1) 



log([OIII]i5007/H8) 



log([NII]^6584/H(!r) - 0.47 

061 

log([NII]i6584/Ha) - 0.05 



+ 1.30 (2) 



log([OIII]/l5007/H/3) = (-30.787 + 1.1358 log([NII]^6584/HQ') 
+0.27297 (log([NII]/l6584/HQr)^))tanh(5.7409 log([NII]i6584/Ha)) 

-31.093 

Figure[T]shows how these boundaries map on to the standard 
emission line diagram for our control sample. We consider galax- 
ies classified by the KOI scheme as AGN, to be 'pure' AGN and 
galaxies classified by S06 as star-forming to be 'pure' star-forming. 
Galaxies classified as AGN by S06 and star-forming by KOI are 
considered as 'composite' (e.g. Ho, Filippenko & Sargent 1993) 
objects as they may contain contributions from both AGN and HII 
region spectra. Galaxies may also reside in the composite zone if 
they have particularly high star-formation rates, high metallicities 
or harbour shock-excited gas. AGN may themselves be divided into 
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Figure 4. The fraction of pair galaxies (solid points) classified as AGN according to the S06 diagnostic as a function of projected separation. Open points 
show the AGN fractions for the control galaxies matched to the pairs in each bin of projected separation. The panels show the full pairs sample (top left), only 
the approximately equal mass (major) pairs (top right), the greater (bottom left) and lower mass (bottom right) in unequal mass (minor) pairs. 



Seyfert galaxies and those that exhibit weaker, low ionization emis- 
sion lines, the so-called 'low ionization nuclear emission line re- 
gion' objects (LINERs, Heckman 1980). Although the excitation 
mechanism for LINERs is itself ambiguous, it is likely to origi- 
nate from a non-stellar process (Ho et al. 1993). In this paper, we 
do not consider these finer AGN classifications or the effects of 
shock excitation. Although the latter may be present in mergers, 
the SDSS spectra can not easily diagnose its presence (Rich et al. 
2010, 2011). 



4 THE AGN FRACTION IN GALAXY PAIRS 

We begin simply, by using the BPT diagram to classify galaxies as 
AGN or non-AGN. For this exercise, we first use the S06 diagnostic 
to distinguish between the two galaxy classes, as it includes even a 
modest amount of nuclear activity. Galaxies with a QSO classifica- 
tion are also included in the AGN sample. Figure|2]shows the AGN 
fraction of pairs as a function of projected separation (r,,). Note that 
we quote the AGN fraction of the full pairs sample, which contains 
both emission line galaxies that can be classified as AGN on the 
BPT diagram, and quiescent galaxies with no strong emission lines. 
For comparison, we show results both with and without the minor- 
ity of galaxies with a QSO classification. There is a steady increase 
towards small projected separation in the AGN fraction for pairs 
with Vp < 40 /j^Q kpc. The control galaxies show no such depen- 



dence on Tp, where the projected separation for the control galaxies 
reflects the value of the pairs to which they are matched. This re- 
sult qualitatively supports the findings of Alonso et al. (2007) who 
found a 10% increase in AGN fraction at r,, < 25 kpc (see also 
Woods & Geller 2007 and Rogers et al. 2009). However, we note 
that the actual fractional increase in AGN is not comparable be- 
tween ours and previous works, whose samples contain only emis- 
sion line galaxies (see Section|6]for further discussion). The lower 
panel of Figure|2]shows the ratio of the AGN fraction in pairs rela- 
tive to the control, i.e. the excess of AGN in the pairs. At the closest 
separations (r^ < 10 /j^q kpc) AGN are 2.5 times more common in 
pairs than in the control sample. The inclusion of QSOs makes a 
negligible impact on these statistics. 

By splitting the sample by redshift at a value of 0.1, we show 
in Figure|3]the AGN excess in two redshift bins. The z > 0. 1 sample 
has a slightly higher AGN excess at all values of r^, but the signifi- 
cance is < Icr except at the smallest separations (r,, < 10 h^^ kpc). 
Any increase in AGN fraction in the higher redshift bin is unlikely 
to be due to the effect of covering fraction for two reasons. First, 
the AGN excess is relative to a mass- and redshift-matched control 
sample, so that we expect the covering fraction to be approximately 
the same between a given pair galaxy and its controls. Second, the 
higher covering fraction of high redshift galaxies should lead to a 
lower fractional contribution from AGN, rather than an enhanced 
AGN fraction. Although Figure |3]may therefore indicate some de- 
pendence of AGN fraction on redshift, the redshift interval probed 
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by our sample is too narrow to allow a rigourous investigation of 
the redshift dependence of the merger- AGN connection. 

In Figure |4] we plot the AGN fraction of paired galaxies as 
a function of projected separation for different mass ratios. Once 
again, the AGN fractions of the control galaxies matched to the 
pairs in each bin of projected separation are shown for comparison. 
The major (equal mass) pairs show an enhancement that is very 
similar to the full sample. However, in minor (unequal mass) pairs, 
the more massive galaxies show a strong increase in the AGN frac- 
tion, whereas the signal is marginal in the less massive companion. 
This mass dependence is the inverse of what has been previously 
observed for triggered star formation in close pairs, where the pri- 
mary galaxy shows no enhancement, but the secondary galaxy does 
(Woods & Geller 2007), as predicted by simulations (Cox et al. 
2008). None the less, our AGN results agree with those of Woods 
& Geller (2007) who report a possible (1-2 a) increased AGN frac- 
tion for the primary galaxy in a minor merger. On the other hand, in 
a study of dual AGN, Liu et al. (201 lb) find that both [OIII] lumi- 
nosities and black hole accretion rates are higher in the less massive 
than the more massive component of a minor pair. 

It is interesting to note that in the secondary galaxy of unequal 
mass pairs in our sample at projected separations r,, > 30 Ji^^ kpc 
there is a hint that the AGN fraction may actually be lower than 
in the mass matched control. The lack of enhanced AGN fraction 
in the secondary galaxies is presumably not due to the lack of fuel 
being supplied to the galactic centre. Simulations predict significnt 
torques in the secondary galaxy (Cox et al. 2008) and SFRs are 
indeed observed to be elevated (Woods & Geller 2007). The relative 
trends of enhanced SFRs and AGN fractions in pairs of different 
mass ratios may therefore provide some clues to the timescales on 
which the two phases dominate. 

In Figure|5]we investigate how the choice of AGN diagnostic 
(KOI, K03 or S06) influences the AGN fraction as a function of r,, 
(this figure necessarily excludes QSOs which do not have BPT clas- 
sifications). As expected, the overall AGN fractions increase from 
KOI to S06. The lower panel of Figure [5] shows the AGN excess 
relative to the control sample. Despite some small bin-to-bin dif- 
ferences, the overall trend for all three diagnostics is an increasing 
AGN fraction towards smaller projected separations with a typical 
over-abundance of AGN at r,, < 10 /jy^ kpc of a factor of about 
two. 

Due to the varying contributions of AGN permitted by the dif- 
ferent diagnostics, the three classifications can be combined to clas- 
sify three populations. 'Pure' AGN and 'pure' star-forming galax- 
ies are based on the KOI and S06 demarcations. Emission line 
galaxies located between the KOI and S06 lines are considered to 
be 'composite' in nature. Figure[6]shows the fraction and excesses 
(relative to the control) of the pure AGN, pure star-forming and 
composite galaxies. All three categories show an excess at small 
separations, indicating that the overall emission line fraction is in- 
creasing in small separation pairs. Presumably the emission line 
population is being fed by relatively quiescent galaxies with suffi- 
cient residual gas that an interaction can spur them into (star for- 
mation or AGN) activity. Interestingly, although one of the best 
known properties of close galaxy pairs is their high star forma- 
tion rates, the actual number of galaxies in the pure star-forming 
category shows the most modest enhancement, increasing by only 
20%. The 'pure' AGN fraction increases by 50%, but the composite 
galaxies, which have significant contributions from both AGN and 
star formation show the largest excess at r,, < 10 kpc, by a fac- 
tor of 2.5. Large fractions of composite (also sometimes referred to 
as low luminosity AGN or transition objects) have been previously 
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Figure 5. The fraction of pair galaxies classified as AGN according to the 
KOI (red), K03 (green) and S06 (blue) diagnostic as a function of projected 
separation. The dashed lines indicate the median AGN fraction in the con- 
trol sample. The top panel shows the absolute fraction and the lower panel 
shows the excess relative to the control. QSOs not present in the MPA/JHU 
line flux catalogues are excluded from this Figure. 



noted in pairs studies by Pastoriza, Donzelli & Donatio (1999) and 
Focardi, Zitelli & Marinoni (2008). 



5 DOUBLE AGN 

In Figure|7]we investigate the fraction of close pairs in which both 
of the galaxies are classified as an AGN. In order to calculate the 
double AGN in a control sample, we make fake pairs by taking the 
best matched control galaxy of each galaxy in a pair. In Figure [7] 
we elect to use the S06 diagnostic, but the same qualitative trend is 
present for the other classifications. Whereas the overall S06 AGN 
fraction for galaxies in pairs (Figure |5j increases by a factor of 
2.5, the fraction of parrs in which both galaxies are AGN is about 
8 times higher than the control. The fraction of double AGN ex- 
pected at random can be estimated from the square of the overall 
AGN fraction in a given r,, bin and is shown with grey points in Fig- 
ure |7] This simple comparison shows that the observed fraction of 
double AGN is higher than expected in the random case by a factor 
of about two and hints at simultaneous triggering in at least some 
of the interacting pairs. Liu et al. (201 lb) have recently used corre- 
lated properties (such as stellar ages, 4000 A break strength and H(5 
absorption equivalent width) to infer synchronous AGN triggering. 

We test for correlated AGN activity within a given pair as a 
function of separation using a statistical approach. For a given frac- 
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Figure 6. The fraction of pair galaxies classified as AGN (red), composite 
(green) and star-forming (blue) galaxies as a function of projected separa- 
tion. The dashed lines indicate the median AGN fraction in the control sam- 
ple. QSOs not present in the MPA/JHU line flux catalogues are excluded 
from this Figure. 



tion of AGN in a population, /, the random (i.e. uncorrelated) frac- 
tions of pairs with zero, one or two AGN (/,<), /ri,/,-2 respectively) 
are: 



=2/(1-/), 

f,-2 = f- 



(4) 
(5) 
(6) 



We now assume that in a given sample with overall AGN (in- 
cluding QSOs) fraction /, some fraction, x, of the pairs have cor- 
related emission line properties, i.e. the occurrence of an identical 
classification is not random but induced by the interaction. There- 
fore, the fraction of uncorrelated pairs is (1 - x). In this scheme, 
the correlated pairs could either be both AGN or both star-forming. 
The fraction of uncorrelated pairs with zero, one or two AGN are 
then obtained by multiplying the random fractions in equations |4] 
-|6]by (1 - x). The fraction of correlated pairs is simply x multi- 
plied by either the AGN fraction, / (for AGN-AGN pairs) or the 
star-forming fraction 1 - / (for SF-SF pairs). Table [T] summarizes 
the equations for correlated and uncorrelated fractions. Note that in 
pairs with one AGN only the uncorrelated case is applicable. Fi- 
nally, the actual observed fraction of pairs with zero, one or two 
AGN will be the combination of the correlated and uncorrelated 
fractions (i.e. the sum of the two columns in Table[TJ: 





CoiTelated fraction 


Uncorrelated fraction 


OAGN 


41-/) 


(l-x)(l-/)2 


1 AGN 




(l-x)2/(l-/) 


2 AGN 







Table 1. Fraction of pairs with zero, one or two AGN that are either cor- 
related or uncorrelated, where / is the overall AGN fraction and x is the 
fraction of correlated pairs. 



/ =(l-x)2/(l-/), 

f2=xf+(l-x)f. 



(8) 
(9) 



The observed fractions of pairs with zero, one or two AGN 
(/o,/i,/2 respectively) and the observed total AGN fraction (/) in 
any given r,, bin can therefore be combined with any of the equa- 
tions|7]-|9]to determine x. For example, from equation|9] 



f2-f 



(10) 



/o = x(l-/) + (l-x)(l-/)2 



(7) 



The fraction of pairs with correlated AGN (x/) is shown as 
a function of r,, in the lower panel of Figure |7] The errors on x/ 
are determined with jackknife re-sampling. In each r,, bin xf is re- 
calculated after systematically removing the AGN pair in turn 
and calculating (5, = xf - x;/. For the A' pairs in each r,, bin the 
error is then [(A^ - l)/A'S,52]i/2. 

As expected, the correlated fraction of the control sample fol- 
lows the zero line. For the pairs, there is an increase in the frac- 
tion of pairs with correlated double AGN as the projected sepa- 
ration decreases. Even pairs with separations of 40 < r,, < 80 
/lyg kpc exhibit a non-zero fraction of double AGN pairs that are 
correlated. The raw AGN fractions in Figure [2] show a strong ex- 
cess of AGN only within 40 hZj^ kpc, with a statistically less robust 
(Icr) enhancement out to 60 /j^q kpc. However, the sample is likely 
to become increasingly contaminated by interlopers (galaxies that 
are not truly interacting) at wider separations. An elevated (but de- 
creasing) AGN fraction at wide separations may therefore become 
increasingly difficult to detect. The results from the double AGN 
indicate that the effect of the interaction on nuclear activity may 
actually be present out to 80 /i^g kpc and beyond. This would indi- 
cate that the AGN phase lasts for a significant fraction of the oribital 
time of the merger. Interestingly, Patton et al. (2011) find that the 
central galaxy colours of pairs remain bluer than their control sam- 
ple out to similarly wide separations. The blue colours can be ex- 
plained by elevated star formation rates in pairs with projected sep- 
arations out to 80 /j^Q kpc (Scudder et al. in preparation). It would 
be of great interest to extend star formation and AGN analyses to 
even wider separations. 



6 DISCUSSION 

6.1 Comparison to previous work 

As reviewed in the Introduction, there has been considerable 
disagreement in the literature regarding the connection between 
galaxy mergers and AGN activity. The two main techniques that 
have previously been used to tackle this connection are 1) through 
the imaging of AGN host galaxies (to search for morphological 
disturbances) and 2) the presence of AGN in merger candidates. 
The former technique has the obvious difficulty of imaging depth 
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Figure 7. Top panel: The fraction of pairs witli two AGN as a function 
of projected separation (filled black points). The grey points indicate the 
expected (random) double AGN fraction based on the values in Figure |2] 
There is a clear excess of observed double AGN relative to the random case. 
A sample of control pairs is constructed by taking the best matched control 
(out of ten) of each paired galaxy and pairing it up with the best matched 
control to the companion galaxy. The fraction of double AGN in the control 
pairs is shown with open symbols. Lower panel: The fraction of correlated 
AGN pairs (xf, see Table[T) as a function of projected separation. Pairs are 
shown with solid symbols, control galaxies are shown with open symbols. 
The enhanced fraction of correlated double AGN (relative to the control) 
indicate that both galaxies in the pair are undergoing synchronized activity. 



and the ability to detect morphological disturbances. These issues 
have been discussed in the literature and we refer the reader to the 
papers focused explicitly on this technique for more details (e.g. 
Ramos Almeida et al. 2011a,b). In this paper, we have taken the 
second approach: to quantify the AGN fraction in a sample of likely 
merger candidates. This technique also has its controversies. Even 
with large samples and carefully matched controls, a disagreement 
remains as to whether close pairs have a higher AGN fraction (Li et 
al. 2006, 2008; Alonso et al. 2007; Woods & Geller 2007; Ellison 
et al. 2008. Darg et al. 2009; Rogers et al. 2009). 

Our results can be used to shed light on some of the reasons 
behind the discrepant results in the literature, as well as some of 
the selection effects which do not seem to play a major role. For 
example, the choice of mass ratio has a relatively minor effect on 
the detection of an enhanced AGN fraction. We find that restricting 
our analysis to approximately equal mass mergers produces quite 
similar results to the full 1:10 mass range (Figure |4]l. The result is 
also robust to the choice of relative velocity: strong enhancements 
in the AGN fraction are still seen at small separations even with 
a relatively relaxed threshold of AV < 500 km s '. The choice 



of AGN diagnostic is also not a major factor; we see increasing 
AGN fractions at smaller projected separations for all three of the 
standard BPT diagnostics (Figure|5}. 

One selection criterion that does play a significant role is the 
S/N imposed on the detection of emission lines. The S/N>5 cut 
used in our work is fairly typical of emission line studies wishing to 
accurately determine gas-phase metallicities (as reviewed by Kew- 
ley & Ellison 2008). However, the location of a galaxy on the BPT 
diagram may be less sensitive to S/N ratio. Moreover, some au- 
thors have argued that imposing S/N cuts biases the sample against 
low luminosity AGN which preferentially have LINER-like spec- 
tra (e.g. Woods & Geller 2007; Cid-Ferandes et al. 2010). We con- 
firm this general trend by calculating the fraction of LINERs in the 
full SDSS sample as a function of S/N. Although LINERs lie in 
a slightly different part of the BPT diagram (e.g. Figure 3 of Cid- 
Femandes et al. 2010) we classify them using the [01] diagnostic 
of Kewley et al. (2006b). The LINER fraction rises from 4 per cent 
for a S/N> 5 to 17 per cent for a S/N>1. However, the fraction of 
AGN in the pairs sample is actually lower for lower S/N cuts. For 
example, at the closest separations (r,, < 10 /i^^ kpc), the AGN ex- 
cess is only 20 per cent for a cut of S/N> 2. When the criterion is 
increased to S/N> 8, the AGN excess is a factor of 3.5 (recall that 
for our fiducial S/N> 5 cut the enhancement is a factor of 2.5). This 
trend of higher AGN fraction for more stringent S/N criteria indi- 
cates that AGN in the pairs sample are not dominated by objects 
with LINER-like spectra. 

One important component of our analysis has been a care- 
ful treatment of selection effects. As discussed in Section |2l in- 
completeness at small separations in the SDSS leads to a biased 
distribution of masses as a function of projected separation. The 
AGN fraction is extremely sensitive to galactic stellar mass: for 
M^, = 10'" Mq the fraction of emission line galaxies classified as 
AGN by S06 is ~ 20%, rising to about 75% by = 3 x 10'° Mq. 
There are similarly steep dependences of the AGN fraction when 
using the KOI and K03 diagnostics, although the overall fraction 
is lower than with S06. Therefore, even a 0. 1 dex bias in mass can 
change the fraction of emission line galaxies classified as AGN by 
10%. We have mitigated this effect through a cull of our wide pairs 
(to remove mass dependences on projected separation) and a care- 
ful matching of a control sample (to facilitate a relative compari- 
son). 

The enhanced AGN fraction at close separations is striking 
when displayed as a function of projected separation (Figure O. 
The significance of our result would be diminished if we considered 
the AGN fraction of the pairs sample as a whole (i.e. not as a func- 
tion of Vp) and compared it to the control. Our work has therefore 
benefitted from a large sample where previous works have calcu- 
lated AGN fractions for pairs and controls in a wholesale way (e.g. 
Alonso et al. 2007; Woods &. Geller 2007; Ellison et al. 2008). It is 
also worth noting that our pairs sample differs from these previous 
works, in that it includes galaxies without strong emission lines. 
The AGN fractions quoted here are fractions of all galaxies, not 
fractions of emission line galaxies that are classified as AGN. Re- 
peating our analysis for an emission line only sample shows a much 
more modest increase in the AGN fraction at small separations, 
by only about 20% (using the S06 diagnostic). As demonstrated 
by Figure [6] the overall fraction of galaxies with emission lines 
(star-forming, composite and AGN) increases in pairs with smaller 
projected separations. Calculating AGN fractions for a complete 
galaxy sample therefore accounts for the migration of quiescent 
galaxies into the emission line regime. Our technique therefore 
tracks a single sample of galaxies and their changing emission line 



AGN in galaxy pairs 9 



properties. For an emission line only sample, a changing AGN frac- 
tion as a function of r,, is then a combination of the changing size 
emission line sample, plus the changing fraction of AGN. 



6.2 Multiple AGN channels? 

The arguments above help to reconcile the disparate results in the 
AGN fraction of close pairs by demonstrating some of the selection 
effects that can lead to different AGN fractions. However, our con- 
clusion that mergers can lead to AGN should be distinguished from 
whether all AGN are the result of mergers. In theory, any process 
that can funnel gas to the centre of a galaxy may be expected to 
increase the accretion rate of the black hole and lead to an AGN. 
Indeed, there is an observational link between enhanced central star 
formation rates and black hole activity that extends beyond mergers 
to the general galaxy population (Cid-Fernandes et al. 2001; Kauff- 
mann et al. 2003a; Silverman et al. 2009; Reichai'd et al. 2009). 
Other evidence has led to claims that mergers are indeed unlikely 
to be the sole AGN trigger (e.g. Li et al. 2006; Lutz et al. 2010; Mul- 
laney et al. 201 1), which may play a role in explaining why many 
imaging surveys of AGN hosts detect no signs of (an excess of) 
recent merger activity (e.g. Grogin et al. 2005; Gabor et al. 2009; 
Cistemas et al 201 1; Kocevski et al. 201 1). Two channels for AGN 
activity are often recognised in simulations, the violent, merger- 
induced 'quasar-mode' and the more gentle 'maintenance-mode' 
(e.g. Di Matteo et al. 2005; Hopkins & Hemquist 2006; Hopkins et 
al. 2008). There is general observational support for the association 
of quasars and radio galaxies with mergers (e.g. Canalizo & Stock- 
ton 2001; Ramos Almeida et al. 2011a,b and references therein) 
when suitable control samples are used. However, the work pre- 
sented here is perhaps more surprising (and challenges the simple 
two channel model described above) in that we find an excess of 
relatively low luminosity (mostly Seyfert-like) AGN in interacting 
pairs. 

If there are multiple pathways to creating an AGN, we should 
distinguish between asking 'what fraction of close pairs have 
AGN?' (the question addressed in the main body of this paper) 
and 'what fraction of AGN are associated with mergers?'. The lat- 
ter question epitomizes the approach of most imaging surveys that 
have not found a significant difference between the disturbed, or 
close pair, fractions in AGN and non-AGN samples. We can also 
ask this second question of our SDSS sample. The same AGN clas- 
sification criteria described in Section[3]are applied to the complete 
SDSS main galaxy catalogue. Of the 43,436 galaxies that are clas- 
sified as AGN, 0.030±0.001 have a close companion within 200 
km s"' and 30 hj^ kpc. We also define a sample of 43,436 non- 
AGN galaxies; the best simultaneous match in stellar mass and red- 
shift to each of the AGN galaxies. The fraction of non-AGN with a 
close companion within 200 km s"' and 30 /i^^ kpc is 0.020±0.001 
(the error is based on simple poisson statisticsjj. Therefore, from 
our sample of pairs we can say that close interactions result in a 
higher AGN fraction and AGN in the full sample show an excess of 
close companions. This contrasts with the imaging studies (mostly 



^ The fraction of galaxies with a close companion in both the AGN and 
non-AGN samples is not accurate in the absolute sense, due to spectroscopic 
incompleteness and other selection biases. For this simple comparison, we 
have not attempted an accurate correction, but assume that both the AGN 
and non-AGN are affected identically and hence the fraction of AGN and 
non-AGN with a close companion can be compared in a relative sense. 



at higher redshift) listed above which do not find an excess of merg- 
ers associated with X-ray selected AGN. However, even accounting 
for the spectroscopic incompleteness (Patton & Atfield 2008), it is 
clear that not all AGN in the SDSS sample have a close companion. 

If mergers are not the dominant cause of AGN activity, then 
the next goal is to identify the main mechanism that leads to black 
hole accretion. Simulations provide some insight here, by investi- 
gating the relative efficiencies of direct 'cold flow' gas accretion 
and mergers (e.g. Hopkins & Hernquist 2006; Di Matteo et al. 
2011). Indeed, Bournaud et al. (2011) argue that, at least at high 
redshift where galactic gas fractions are high due to inflows from 
cold streams, disk instabilities could lead to significant AGN fu- 
elling. This argument is supported by the prevalence of clumpy 
disks in high redshift galaxies (Elmegreen et al. 2007). However, 
Wisnioski et al. (2011) do not find any AGN dominated disks in 
their small sample of clumpy z ~ 1 disks. Alternatively, Ellison et 
al. (201 1) have recently demonstrated that barred galaxies with log 
M<- > 10'" Mq have SFRs that are enhanced at a similar level as 
close pairs. However, due to their relative ubiquity, the enhanced 
SFRs in barred galaxies contribute at least three times more to the 
central stellar mass build-up than mergers. We will investigate the 
AGN fractions in barred galaxies in a forthcoming paper (Nair et 
al. in preparation). 



7 CONCLUSIONS 

We have compiled a large (11,060) sample of close galaxy pairs, 
selected with careful attention to the photometry at small separa- 
tions, emission line quality and biases due to small separation in- 
completeness. The AGN fraction of the pairs sample is compared 
to a mass- and redshift-matched control sample with 10 controls 
per pair galaxy. 

The main conclusions of this paper are: 

(i) The AGN fraction in close pairs of galaxies increases by a 
factor of up to ~ 2.5 at projected separations r,, < 10 Ii^q kpc. 
Although the increase depends slightly on the choice of AGN di- 
agnostic, elevated AGN fractions are found for closely separated 
pairs independent of the diagnostic choice. 

(ii) The excess of AGN is larger for more stringent S/N cuts, 
indicating that the AGN in pairs are not dominated by objects with 
LINER-like spectra. 

(iii) We discuss how AGN are likely to be triggered by multi- 
ple processes, of which mergers are just one channel. However, in 
contrast to models, our results show that mergers can result in an 
increase in relatively low luminosity AGN, as well as the power- 
ful radio galaxies and quasars that have been investigated by other 
studies. 

(iv) The fraction of AGN in close pairs depends mildly on the 
ratio of stellar masses. The highest enhancements in AGN fraction 
are seen for the equal mass pairs. The lower mass galaxies in un- 
equal mass pairs show relatively little enhancement AGN fraction. 

(v) In addition to an increase in the 'pure' AGN fraction at de- 
creasing projected separations, both the star-forming and composite 
fractions also increase. Increases in all three emission line classes 
indicates that some quiescent galaxies have been transformed dur- 
ing the interaction process. 

(vi) The fraction of pairs in which both galaxies are classified 
as AGN by the Stasinska et al. (2006) AGN diagnostic is a factor 
of up to 8 higher than a sample of 'control pairs'. This exceeds 
the expected random fraction of double AGN in pairs by around a 
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factor of two, indicative of correlated AGN triggering between the 
companions. 

(vii) The correlated double AGN fraction is investigated using a 
statistical approach and found to increase at small separations, but 
remains elevated out to 80 h:^^ kpc. This is evidence that the AGN 
phenomenon may persist in relatively wide separation pairs. 

The results of this paper are in clear support of interaction- 
driven AGN activity that occurs well before final galaxy-galaxy co- 
alescence. Since the merger timescale may be as long as a gigayear 
(Kitzbichler & White 2008), this does not provide a strong con- 
straint on a delay between triggered star formation and AGN activ- 
ity. However, the shape of the double AGN fraction as a function of 
projected separation (Figure|7]l is very similar to the colour offsets 
(Patton et al. 201 1) and star formation rate enhancements (Scudder 
et al. in preparation) in our sample. Scudder et al. compare the ob- 
served star formation rates with simulations and demonstrate that 
they can be reproduced with a burst of star formation occurring 
after the first pericentric passage, followed by a more significant 
burst at final coalescence. The similarity between star formation 
and AGN enhancements as a function of projected separation may 
be indicative of a similar process for accretion onto the nucleus with 
a relatively small time lag. Delays of a few hundred million years 
(Wild et al. 2010; Hopkins 2011) between peaks in star formation 
and accretion rate correspond to changes in projected separation 
of order the bin width in our projected separation figures (depend- 
ing on relative velocity and the orientation). Megayear delays are 
therefore relatively prompt compared with the resolution provided 
by the observations presented here. The high fraction of composite 
galaxies also demonstrates that the AGN and star formation pro- 
cesses are on-going concurrently, with no need to shut down the 
latter before the former begins (c.f. Li et al. 2008; Schawinski et al. 
2009). 
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